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Abstract 

We perform semi-classical molecular dynamics simulations of screening by bound 
electrons in low energy nuclear reactions. In our simulations quantum effects cor- 
responding to the Pauli and Heisenberg principle are enforced by constraints. In 
addition to the well known adiabatic and sudden limits, we propose a new "dissi- 
pative limit" which is expected to be important not only at high energies but in 
the extremely low energy region. The dissipative limit is associated with the chaotic 
behavior of the electronic motion. It affects also the magnitude of the enhancement 
factor. We discuss also numerical experiments using polarized targets. The derived 
enhancement factors in our simulation are in agreement with those extracted within 
the i?-matrix approach. 
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Dynamics; Classical Chaos 



1 Introduction 



The relation between the tunneling process and dynamical chaos has been 
discussed with great interests in the field of nonlinear science Though 
tunneling is a completely quantum mechanical phenomenon, it might be in- 
fluenced by classical chaos. In the sense that chaos induces fluctuations on 
the classical action which essentially determines the tunneling probability. We 
study the phenomenon by examining the screening effects by bound electrons 
in low energy fusion reactions, where the experimental cross sections with gas 
targets show an increasing enhancement with decreasing bombarding energy 
with respect to the values obtained by extrapolating from the data at high 
energies (0; 0) ■ Many studies attempted to attribute the enhancement of the 
reaction rate to the screening effects by bound target electrons (jl). In this 
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context one often estimates the screening potential({7 e ) as a constant decrease 
of the barrier height in the tunneling region through a fit to the data. A puzzle 
has been that the screening potential obtained by this procedure exceeds the 
value of the adiabatic limit, which is given by the difference of the binding 
energies of the united atoms and of the target atom and it is theoretically 
thought to provide the maximum screening potential (0). The experiments 
have been performed using deuterated metallic targets, as well, and it is re- 
ported that one observes systematically large screening potentials for various 
kinds of metals (M). 



The difficulty lies in the determination of the absolute value of the bare cross 
section. The extrapolated bare cross section, i.e., the astrophysical S-factor, 
might be affected by the choice of the high energy region where one assumes 
the screening potential can be neglected (j8(). Over these years, several ways to 
determine the bare cross sections have been proposed. There are theoretical 
attempts, using the i?-matrix theory (0), and experimental, using the Trojan 
Horse Method(THM) Jl(J[ll|). However, the THM gives only relative values 
of the S-factors and one needs to normalize to direct methods data obtained 
at high incident energies. Thus the comparison between newly obtained bare 
cross sections and the cross sections by the direct measurements gives a variety 
of values for the screening potential. These values are often smaller than the 
sudden limit or larger than the adiabatic limit. Theoretical studies performed 
using the time-dependent Hartree-Fock(TDHF) scheme (flil ; [kJ) suggest that 
the screening potential is between the sudden and the adiabatic limits. 



One of the aims of this paper is to try to assess the effect of the screening 
quantitatively. Up to now, the dynamical effects of bound electrons have been 
studied only in some limited cases with a few bound electrons (the D+d with 
atomic target (fl2[ fl3h and molecular D 2 target (14), the 3 He+(i (fl3 )) with 
the TDHF method. We investigate here the dynamical effects, including the 
tunneling region, for other systems with many bound electrons; D+D, 3 He+D, 
looking the effect of the electron capture of projectile, 6 Li+cZ, 6 Li+D. 

To simulate the effects of many electrons, we use the constrained molecular 
dynamics (CoMD) model Q, Q M). At very low energies fluctuations are 
anticipated to play a substantial role. Molecular dynamics contains all pos- 
sible correlations and fluctuations due to the initial conditions(events). The 
prescription using constraints for the Heisenberg uncertainty principle and the 
Pauli exclusion principle is based on the Lagrange multiplier method (0). In 
extending the study to lower incident energies, we would like to stress the con- 
nection between the motion of bound electrons and chaos. In fact, depending 
on the dynamics, the behavior of the electron(s) is unstable and influences 
the relative motion of the projectile and the target. This feature is caused 
by the nonintegrability of the iV-body system(A r > 3) and it is well known 
that the tunneling probability can be modified by the existence of a chaotic 
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environment. In this connection we propose a new "dissipative limit" which 
is driven by electron ejection. We examine also numerical experiments using 
polarized targets. 



The paper is organized as follows. In sect. 2 we introduce the enhancement 
factor f e and describe the essence of the constrained molecular dynamics ap- 
proach briefly. In sect. 3 we mention limiting cases where one can easily esti- 
mate the screening potential. We discuss particular reactions, the effect of the 
electron capture by projectile and the numerical experiments using polarized 
targets. We summarize the paper in sect. 4. 



2 Formalism 



We denote the reaction cross section at incident energy in the center of mass 
E by cr(E) and the cross section obtained in absence of electrons by cr (E). 
The enhancement factor f e is defined as 



f _ <E) 

Je 



If the effect of the electrons is well represented by the constant shift U e of the 
potential barrier, o~(E) is replaced by a (E + U e ) following |l2l 17), 



fe ~ exp 



where r/(E) is the Sommerfeld parameter (0). We calculate the enhancement 
factor through the CoMD simulation. 



2. 1 Constrained Molecular Dynamics 



We describe the essence of the CoMD briefly. The total Hamiltonian is written 
down as 

N 

if (r; p) = + U(Ti) - m lC 2 ), (3) 
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where we use relativistic kinematics: £j = ypfc 2 + m 2 c A , rrii and 

U{Ti)= £ (4) 
j(#0=o |r " J r ^ 



are the energy, the mass and the potential of the i-ih particle, i.e., an electron 
or a nucleus, respectively. Here, g; is the charge of the particle i. In the classical 
molecular dynamics(CMD) one solves the Hamilton equations, i.e.,: 



dvi pic 2 dpi 

The initial configurations of the bound electrons are prepared using the CoMD, 
as well |l6f). To take the feature of the Pauli blocking into account, we use the 
Lagrange multiplier method for constraints. Our constraints which correspond 
to the Pauli blocking is fa < 1 in terms of the occupation probability and can 
be directly related to the distance of two particles, i.e., r^p^, in the phase 
space. Here = |rj — and = |p* — Pj|. The relation /j < 1 is fulfilled, if 
fijPij > CpfoSsusj, where £p = 27r(3/47r) 2 / 3 2 1//3 , i,j refer only to electrons and 
Si,Sj(= ±1/2) are their spin projection. We can easily extend the approach 
to the Heisenberg principle where the constraint is expressed as r^pij > 
and = 1, i and j refer to the electrons and to the nuclei. Using these 
constraints, the Lagrangian of the system can be written as 



where Af and Af are Lagrange multipliers for Pauli and Heisenberg principles 
respectively. The variational calculus leads to: 



dTj pjc 2 1 Af Af \ dp^ 

dt + hf$[tH 6> Si ' Sj ) 11 [0 

rf P* V7 TTt \ 

^ = -v rf /( rj) - 

where Af and Af are zero, if the two particles i and j are far separated; 
r ijPij > £^ (£ = £p f° r electrons i and j, which have identical spin projections 
and £ = for the other pairs of particles). In this case the eqs. (7) and (8) 
redeem the eqs. (5). The integration of the eqs. (7) and (8) is performed using 
Hermite integration scheme which is efficient and capable of high precision. 
The scheme adopts variable and individual time-steps for each electron and 




(6) 
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nucleus (19). In this way we obtain many initial conditions which occupy 
different points in the phase space microscopically. The Lagrange multipliers 
are used to obtain the ground states of atomic nuclei. During the collision we 
set the Aj = to conserve energy. If needed, a collision term among electrons 
is introduced to avoid the over occupancy of the phase space (|15T). For the 
collisions reported here we found the violation of the Heisenberg and Pauli 
principles negligible during the collision dynamics. 

The importance of the influence of the tunneling region to the electron screen- 
ing has been discussed in fllih . In order to treat the tunneling process in the 
CoMD framework, we define the collective coordinates R coW and the collective 
momentum p co11 as 

R coll = rp _ Tt . pcoll _ pp _ (g) 



where r T ,r P (pr, Pp) are the coordinates(momenta) of the target and pro- 
jectile nuclei, respectively. When the collective momentum becomes zero, we 
switch on the collective force, which is determined by F c p U = p co11 and = 
— P coH , to enter into imaginary time (j2oh. We follow the time evolution in the 
tunneling region using the equations, 



dr T(P) _ Pt(P) . d PT(P) _ ^ , :,„0s , ,,r-H 

dr St(p) ' dr 



V r U(4 {P) )-2F^, (10) 



where r is used for imaginary time to be distinguished from real time t. r®, p s 
and p%p) are position and momentum of the target (the projectile) during the 
tunneling process respectively. In fact in the tunneling region the path which 
gives an important contribution to the action of Feynman path integral lies in 
the imaginary time region. The path coincides with the classical path in the 



potential valley which is given by turning the barrier upside down (|2lt). 
The penetrability of the barrier is given by (J20I ) 

U(E) = (l + exp(2A(E)/h)y 1 , (11) 

where the action integral A(E) is 

A(E) = J P coU dR coU , (12) 



r a and are the classical turning points. The internal classical turning point 
r b is determined using the sum of the radii of the target and projectile nuclei. 
Similarly from the simulation without electron, we obtain the penetrability 
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of the bare Coulomb barrier Hq(E). Since nuclear reaction occurs with small 
impact parameters on the atomic scale, we consider only head on collisions. 
The enhancement factor is thus given by eq. (1), 



for each event in our simulation. Thus we have an ensemble of f e values at 
each incident energy. 



3 Applications to Electron Screening Problem 

We investigate the enhancement factor for the reactions D+d, D+D, 3 He+<i, 
3 He+D, 6 Li+d, 6 Li+D. It is well known that in the low incident energy region 
the projectile often captures electrons before it collides with the target nucleus. 
The atomic projectile cases(D+D, 3 He+D, 6 Li+D) in contrast with the bare 
ionic beam projectile cases(D+d, 3 He+d, 6 Li+d) reveal the effect of the extra 
electrons. 

3.1 Two well known limits and the Dissipative limit 

As it is often discussed in the literature, there are two limiting cases where one 
can easily estimate the screening potential U e . One is the adiabatic and the 
other is the sudden limit (jD) . In the sudden limit one assumes that the relative 
velocity of the two ions is relatively faster than that of the electrons, so that the 
bound electrons remain frozen during the reaction. Instead in the adiabatic 
limit which is associated with the slow relative velocity limit, the electrons 
change their configurations following the relative motion of ions. Especially 
in the case where the electrons occupy the ground state of the system at the 
beginning of the reaction, the electrons tend to continue occupying the ground 
state throughout the reaction process. Effectively the electrons move close to 
the target nucleus during the collision. In such a way the electrons continue 
to shield the coulomb repulsive field between the target and the projectile. 
In this limit the screening potential is given by the difference between the 
binding energies of the target (B.ES T ^) and the united atom (B.ES UA ') which 
is formed in the last stage of the reaction, 



We propose another case which is obtained as a result of the bound electron 
emission. It gives usually a negative effect to the fusion enhancement, because 



f e = U(E)/U (E) 



(13) 



U( A °) = B.E.W - B.ES UA \ 



(14) 
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the electrons take the kinetic energy away from the ions relative motion. We 
name the new limit as "dissipative limit" (DL). The screening potential in the 
dissipative limit is given by 

= B.ES T \ (15) 

where we assume that the ejected electron has zero kinetic energy. If the elec- 
tron takes finite kinetic energy away, the screening potential can be smaller 
than the dissipative limit. Note that U^ DV> has a negative value. Such a phe- 
nomenon seldom occur in the framework of the mean field approach such as 
TDHF because of the usually assumed spherical (or cylindrical) symmetry of 
the system. We would like to stress here that the DL could be used to per- 
form low energy collisions using relatively higher energy beams, similar in some 
sense to the THM. In fact we propose to perform experiments where the fusion 
residues are detected in coincidence with energetic electrons. Thus increasing 
the value of the ejected kinetic energy of the electrons is as if decreasing the 
beam energy. 

3.2 Reactions 

3.2.1 D+d and D+D reactions 

Fig. 1 shows the incident energy dependence of the enhancement factor for 
the reactions D+d and D+D, where the systems involve 1 and 2 electrons 
respectively. The open and closed squares show the average enhancement fac- 
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Fig. 1. Enhancement factor as a function of incident center-of-mass energy for the 
reactions D+d and D+D. Error-bars represent the variances obtained from the 
events generated for each beam energy. 

tors f e over events for the reactions D+d and D+D, respectively. The vari- 
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ances £ = y/ e 2 — (/ e ) 2 are shown with error bars. The dotted and dash- 
dotted curves show the enhancement factors in the adiabatic limit fj^ AD ^ for 
an atomic deuterium target and it is obtained by assuming equally weighted 
linear combination of the lowest-energy gerade and ungerade wave function 
for the electron, reflecting the symmetry in the D+d, i.e., 

f(AD) = 1 L^(E)\- + e ^B)\- j ) (16) 



where = 40.7 eV and U^> = 0.0 eV (Q Q for D+d case. If we take 
into account the electron capture of the projectile, i.e., in the case of D+D, 
the enhancement factor in the adiabatic limit is 

1 ,,(s-s.) q r/ (les) 



where [/QmO = 51.7 eV and £/ e (les) = 31.9 eV 0). The solid curve and dashed 
curve show the enhancement factors in the dissipative limit f^ DL ^ for the re- 
actions D+d and D+D respectively. Notice how the calculated enhancement 
factor with their variances nicely ends up between the adiabatic and the dis- 
sipative limits. We performed also a fit of our data using eq. (2) including the 
very low energy region and obtained U e = 15.9 ± 2.0 eV for D+d case and 
U e =21.6 ± 0.3 eV for D+D. 

In Fig.l we saw the rough property of the enhancement factors as aver- 
age values and their variances. Now we discuss their distributions for the 
reaction D+d. Fig. 2 shows the histograms of the normalized enhancement 
factors(/ e // e ) at the incident energies E cm =0.25keV(top panel), 2.5keV(middle 
panel) and lOOkeV (bottom panel). The stars, open squares and crosses specify 
the adiabatic, sudden and dissipative limits at each incident energy, respec- 
tively. A remarkable feature in those figures is that at the highest incident 
energy, the enhancement factor for each event in an ensemble distribute al- 
most as a ^-function around the average enhancement factor. As the incident 
energy goes down, it spreads around the average and the adiabatic limit at 
E cm =2.5keV. However, at the lowest incident energy E cm =0.25keV, where 
any experiment has not reached yet, instead one can see clearly that the dis- 
tribution of the enhancement factor changes and the peak of the distribution 
is in between the sudden and the dissipative limits. At a glance it appears 
strange that the enhancement factor goes to the sudden and dissipative lim- 
its (which we expected to be important in the high incident energy region). 
The phenomenon is due to the chaotic behavior of the bound electrons in 
the many-body systems (0; H^), which becomes important especially in the ex- 
tremely low energy reactions, where the interaction time is relatively long and 
the electrons have enough time to interact with the projectile ions. 
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Fig. 2. Histograms of the normalized enhancement factor for the reaction D+d at 
various incident energies. In the figure stars, open squares and crosses refer to the 
adiabatic, sudden and dissipative limits, respectively. 

To strengthen the finding of the electrons influence on the enhancement factor, 
in Fig. 3 we show the enhancement factors for many events as a function of the 
electrons binding energy at the external classical turning point. As one can 
see, there is a clear correlation. In fact the larger the electrons binding energy, 
the higher the enhancement factor is. Thus in order to control the fusion 
probability we need to have some control on the electronic motion especially 
at low energies as we describe next. 



3.2.2 Polarized targets 

We prepared numerically ensembles of target atoms where the electrons motion 
is polarized perpendicular or parallel Pn to the beam axis. Fig. 4 shows 
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Fig. 3. The enhancement factor as a function of the binding energies of the elec- 
tron at the external classical turning point divided by the incident energies for the 
reaction D+d 

the initial configuration dependence of the enhancement factor at the incident 
energy E cm =1.25keV. The open circles are the enhancement factors for 
targets and the closed squares are the ones for Py targets. It is clear that 
the P± targets give always large enhancement factors and a small variance, 
as it is shown with the error bars in the figure. Instead the Pii targets give 
relatively small enhancement factors and large variances. As we mentioned in 
the reference (jU) or it is discussed in the reference (jj^h as an example of the 
coplanar collisions, this is the case where the oscillational force mostly affects 
the relative motion between the two nuclei. Because of the non-integrability 
of the system under such a force, the motion of the electron becomes unstable 
and it is often ejected to the continuum state. Notice that it corresponds to 
the dissipative limit which we denned in this paper. The large variances of 
the P|| targets originates from the fact that the chaotic behavior of the system 
affects the determination of the enhancement factor. A remarkable feature is 
that with the Py targets the enhancement factor often becomes less than 1. It 
means that in this case the bound electron gives an hindrance to the tunneling 
probability. 

In Fig. 5 we show the incident energy dependence of the average enhancement 
factor for the Pj_ and Py targets with pluses and crosses, particularly in the 
low energy region. The enhancement factors from the P^ targets are always 
larger than that from the Py targets. In contrast to the average enhancement 
from the Pj_ targets, which increases monotonically as the incident energy 
becomes smaller, the average enhancement from the Py targets fluctuates. It 
has also large variances at low energies. 

In figures 6 and 7 we show the histograms of the enhancement factor for Py (top 
panel), P^ (bottom panel) at two relatively low incident energies(.E cm = IkeV 
and 0.25keV). In the abscissa the enhancement factor is normalized with its 
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Fig. 4. The initial configuration dependence of the enhancement factor at the inci- 
dent energy E cm =1.25keV. The enhancement factors for Pj_ targets(open circles) 
and for Pm targets (closed squares) are shown. 
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Fig. 5. same as Fig. 1 but for Pj_ and Pm targets. 



average value over the ensemble of unpolarized targets. 

One notices that at the incident energy E cm = IkeV (Fig. 6) the enhancement 
factor for each event distributes around and close to its average value, which 
is slightly less than the adiabatic limit, as a comprehensive feature. Looking 
carefully at each panel, one sees that the Py target gives smaller enhancement 
factors than the P^ target, as it is shown in Fig. 4. There is also some hindrance 
in the parallel target case. 



The change from Fig. 6 to Fig. 7 is drastic. First, all the events are smaller than 
the adiabatic limit, as we already saw in Fig. 2 for the case of the unpolarized 
targets. Also using P^ targets we obtain enhancement factors always larger 
than the average value. 
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Fig. 6. Histograms of the normalized enhancement factors for the reaction D+d. 
Incident energy E cm =l.keV, polarized parallel target to the beam axis(top panel) 
and polarized perpendicular target (bottom panel). Again stars, open squares and 
crosses indicate the adiabatic, sudden and dissipative limits, respectively. 
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Fig. 7. same with Fig. 6 but at the incident energy E cm =0.25keV. 
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3.2.3 3 He+d and 3 He+D reactions 



An excess of the screening potential was reported for the reactions 3 He+<i with 
atomic gas 3 He target, and D 2 + 3 He with deuterium molecular gas target, for 
the first time in (J3J). Since then various experiments have been performed for 
these reactions. The incident energy covers from 5 keV to 50 keV for 3 He+cL 
Though once the problem of the discrepancy between experimental data and 
theoretical prediction seemed to be solved by considering the correct energy 



loss data (J2J), recent measurements using measured energy loss data (|2f 



report larger screening potentials than in the adiabatic limit for both reactions. 

The electron capture by the projectile plays a minor role in the case of 3 He+d, 
since electrons are more bound in helium targets. However in the recent mea- 
surement Aliotta et al. was performed using molecular Dg" and Dg" targets (25). 
Thus we assess the contribution from the reaction 3 He+D, as well. 

The enhancement factor in the adiabatic limit give f/ e =119 eV for 3 He+<i and 
£/ e =110 eV for 3 He+D, respectively. These are shown in the figure 8 with the 
solid curve for 3 He+d and with the dashed curve for 3 He+D. The comparison of 
these two adiabatic limits implies that the electron capture of projectile would 
give a hindrance compared with the bare deuteron projectile. Meanwhile the 
latest analysis of the experimental data using i?-matrix two level fit (9) sug- 
gests the screening potential U e = 60 eV(corresponding enhancement factor is 
shown with dotted curve). The comparison between direct measurement and 
an indirect method, the Trojan Horse method, suggests the screening poten- 
tial U e = 180±40 e V ( the corresponding enhancement factor is shown with 
dot-dashed curve) (26). The average enhancement factors f e over events in 
our simulations using the CoMD are shown with the open and closed squares 




1 10 
Ecm[keV] 



Fig. 8. Enhancement factor as a function of incident center-of-mass energy for the 
reactions 3 He+(i and 3 He+D. 

for the reactions 3 He+<i and 3 He+D, respectively. The enhancement factors 
of the both reactions 3 He+<i and 3 He+D are in agreement with the extracted 
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values using the i?-matrix approach within the variances over all the events. 
Notice that our calculated enhancement factors for the two systems display 
an opposite trend as compared to the adiabatic limits. The average enhance- 
ment factor of the reaction 3 He+D agrees with the estimation of the adiabatic 
limit and the reaction 3 He+d is below the corresponding adiabatic limit. The 
paradoxical feature comes from the fact that an electron between the two 
ions is often kicked out during the reaction process, i.e., the electron config- 
uration seldom settles down the 5 Li + ground state in the reaction 3 He+d. It 
is known as autoionization in the context of the Classical Trajectory Monte 
Carlo method (|27l). Instead in the case of the 3 He+D, the deuterium projec- 
tile brings its bound electron in a tight bound state around the unified nuclei 
of 3 He and d, practically it ends up with a ground state configuration of the 
5 Li atom. The fits of the obtained enhancement factors suggests the screening 
potentials U e = 82.4 ± 1.9 eV for the 3 He+d and U e = 102.8 ± 3.0 eV for the 
3 He+D. 



3.2.4 6 Li+d reaction 

The S-factors for the reactions 6 Li+rf, 6 Li+p and 7 Li+p were measured over 
the energy range 10 keV < E cm < 500 keV by Engstler,et al. They used 
LiF solid targets and deuteron projectiles as well as deuterium molecular gas 
targets and Li projectiles. 

In the case of LiF target which is a large band gap insulator, one often approx- 
imates the electronic structure of the target 6 Li( 7 Li) state by the 6 Li + ( 7 Li + ) 
with only two innermost electrons. Thus for all three reactions one expects 
the screening; potential in the adiabatic limit U^ AD) = 371.8 - 198.2 ~ 174 
eV. Instead if one uses the ground state of the 6 Li( 7 Li) atom and of the bare 
deuteron target as the initial state, U^ AD ^ =186 eV (0), which is given by the 
solid curve in Fig. 9 . 

However one should be aware that the deuteron or hydrogen projectile plau- 
sibly moves with a bound electron in LiF solid insulator target (j28h. Un- 
der such an assumption we could estimate the screening potential XJ[ AD) = 
389.9 — 198.2 ~ 192 eV. In the case of molecular D 2 or H 2 gas targets, as well, 
we should consider the electron capture by the lithium projectile. 

The bare S-factors for the same reaction have been extracted using an indirect 
method, the Trojan-Horse Method through the reaction 6 Li( 6 Li,aa) 4 He (11). 



The comparison between direct and the indirect methods gives the screening 
potential U e = 320±50 eV. The corresponding enhancement factors are shown 
with the dash-dotted curve. The contrast between the direct measurement 
data and the theoretical estimation for the bare S-factor using the i?-matrix 
theory gives f/ e =240 eV. It is shown with dotted line. The extracted U e with 
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the two different methods are larger than the adiabatic limit. 



We simulate Li+d and Li+D cases only. In the figure 9 the open and closed 
squares show the enhancement factor for the reactions 6 Li+d and 6 Li+D, 
respectively. 




Fig. 9. same as Fig. 8 but for the reactions 6 Li+d and 6 Li+D. 

Again the average enhancement factors of the reaction 6 Li+D are larger than 
those of the 6 Li+<i. The enhancement factors of the reaction 6 Li+D are in 
agreement with the extracted values using the i?-matrix approach within the 
variances over all the events. The fit of the obtained average enhancement 
factors suggests the screening potentials U e = 152.0 ± 9.9 eV for 6 Li+d and 
U e = 214.4±18.5 for 6 Li+D. The screening potential for the reaction 6 Li+d in 
our simulation does not exceed the adiabatic limit nor extracted values using 
the i?-matrix theory and THM, but one for 6 Li+D verges on the extracted 
values using the .R-matrix approach. 

In Table 1 we summarize our results of the screening potentials for each re- 
action and compare them with the extracted value using the i?-matrix theory 
and THM. The screening potentials are obtained by fitting average enhance- 
ment factors for each reactions and we neglect its variances over all events in 
this fitting procedure. 



4 Summary 



We discussed the effect of the screening by the electrons in nuclear reactions 
at the astrophysical energies. We performed molecular dynamics simulations 
with constraints and imaginary time for the reactions D+d, D+D, 3 He+<i, 
3 He+D, 6 Li+d, 6 Li+D. For all the reactions it is shown that both the aver- 
age enhancement factors and their variances increase as the incident energy 
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Table 1 

Comparison between the screening potential(C/ e ) in the adiabatic limit and in our 
CoMD simulation. The errors in the first column originate from the fitting proce- 
dure. 



^e(keV) 

present results Ref. (9)* THM adiabatic limit 



B+d 15. 9± 2.0 8.7, 7.3 22.0* 

D+D 21.6± 0.3 37.1* 

3 He+d 82.4± 1.9 34, 60, 200 180±40 119 

3 He+D 102.8± 3.0 110 

6 Li+cZ 152.0 ± 9.9 259, 248 320±50 186 

6 Li+D 214.4±18.5 



* the column of Ref. (-9) contain not only the results obtained within the the R- 
matrix approach but also with a polynomial fitting. 

* Ref. J3) 

becomes lower. Using bare projectiles we obtained the average screening poten- 
tial smaller than the value in the adiabatic limit for all reactions. It is because 
of the excitation or emission of several bound electrons during the reactions. 
The comparison between bare and atomic projectile cases for each reactions 
revealed that the electron capture of the projectile guides to larger enhance- 
ments. The derived enhancement factors in our simulation are in agreement 
with those extracted within the i?-matrix approach including the variances 
over all the events. 

We performed numerical experiments using polarized targets for the reaction 
D+d. Using Pj_ targets we obtained relatively large enhancements with small 
variances, instead Pii target gives large variances of the enhancement factors 
and relatively small averaged enhancement factors. It is because with the Py 
targets the force exerted from the electron to the relative motion of the nuclei 
is oscillational, in the direction of the beam axis, and the motion of the electron 
becomes often excited or unstable. It is the case where the chaoticity of the 
electron motion affects the tunneling probability and at the same time the 
enhancement factor of the cross section. 
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